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The promising thermoelectric material TiS2 can be easily chemically doped and intercalated. 
We present here studies of single crystals that are intercalated with excess Ti or Co, or substi¬ 
tuted with Ta. We demonstrate the intrinsic impact of these dopants on the thermal transport 
in the absence of grain boundary scattering. We show that Ta doping has the greatest impact 
on the thermal scattering rate per ion added, leading to a five-fold reduction in the lattice 
thermal conductivity as compared to stoichiometric single crystals. 


I. INTRODUCTION 

Due to the rich variety of their physical properties, 2D 
materials are attracting much attention. Huge thermal 
conductivity values (> 10^ W in graphene^ and 

in boron nitride^, various electronic properties in tran¬ 
sition metal dichalcogenides^, high-Tc superconductiv¬ 
ity in cuprates^ are among some of the specific proper¬ 
ties of interest. About two-thirds of the MX 2 transition 
metal dichalcogenides^ crystallize in layered structures 
and they exhibit a wide range of tunable electronic states 
(semiconducting, insulating, semi-metallic and supercon¬ 
ducting), as well as exotic optical and mechanical prop¬ 
erties. Their crystals can be easily cleaved as successive 
X layers are only weakly bonded by van der Waals forces. 

The well-known sulfides derived from TiS 2 have re¬ 
cently been under active investigation in the new con¬ 
text of thermoelectric materials research^^— . While the 
thermoelectric power factor, Pp = S'^/p, {p and S are the 
electrical resistivity and Seebeck coefficient, respectively) 
of TiS 2 single crystals is similar to that of Bi 2 Te 3 , the 
thermoelectric material most commonly used in appli¬ 
cations, the thermal conductivity (k) exhibits relatively 
high values^: K(300Ar) = 6.8WK“^m“^, compared to 
1.2WK“^m“^ for Bi 2 Te 3 iS. This reduces the dimen¬ 
sionless figure of merit, ZT = PpT/n, which measures 
the theoretical efficiency of energy conversion. Several 
attempts have been made to reduce k in dense ceramics, 
for example via intercalation in the van der Waals gap as 
in Cu 3 ,TiS 2 , or substitution at either anion or cation sites 
as in TiS 2 -ySey or Ta 2 Tii_;;S 2 ^“— . These investigations 
have lead to ZT values of 0.4 at 650 K. 

However, for all these ceramic series, it is difficult to 
distinguish whether changes in k are a result of changes 
in the charge carrier concentration, the concentration 
of substituted or intercalated ions, or grain boundary 
effects. Moreover, the chemical formula should more 
properly be written as non-stoichiometric Tii+a;S 2 . It is 
known to be a degenerate extrinsic semiconductor with 
extra Ti^+ cations in the van der Waals gap, the donated 
electrons of which populate the empty Ti d-band of the 


TiS 2 layer—. Meanwhile, the crystallographic textura- 
tion of ceramic samples implies that the transport prop¬ 
erties always result from a convolution of in-plane and 
out-of-plane contributions, not to mention scattering by 
grain boundaries. 

In this article, we report on in-plane electronic 
and thermal transport for several metallic single crys¬ 
tals derived from TiS 2 . The changes in k (in-plane) 
are dominated by the lattice part with a peak at 
0 _d/ 5 (0D=Debye Temperature) expected for phonon- 
dominated heat transport. We show that reduction in the 
lattice heat transport [Kiatuce) deriving from the substi¬ 
tution or intercalation of foreign cations accounts for the 
majority of the change in k, and we identify the species 
that are most effective at reducing the thermal conduc¬ 
tivity. The Kiattice(300K) values obtained for crystals 
derived from TiS 2 can be reduced by a factor of around 
five, from 6.3WK“^m“^ (for the most stoichiometric 
crystals^) to 1.2WK“^m“^ (this study). We conclude 
that the thermoelectric properties of the 2D dichalco- 
genides could be improved by thermal conductivity engi¬ 
neering, even before nanostructuring is pursued. 


II. METHODS 

Large platelet-like crystals were grown by a vapour 
transport technique using iodine as the transport 
agenti^. Intercalated and chemically substituted single 
crystals were grown by pre-reacting Co and Ta with TiS 2 
precursor at 500 °C prior to the growth phase. Shiny thin 
crystals were obtained for all the different batches and 
powder X-ray diffraction confirmed the presence of the 
TiS 2 phase. Scanning electron microscopy with coupled 
EDX chemical analysis on all samples confirmed the nom¬ 
inal compositions within an uncertainty of ±5%. In the 
case of Coo.o 4 TiS 2 , this means that Co was detected but 
with greater error than the nominal composition. One ce¬ 
ramic sample of TiS 2 grown with a slight nominal excess 
of Ti and densified by Spark Plasma Sintering was also 
used for comparison to the single crystals. Full details of 
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FIG. 1. a) Electrical resistivity of the TiS 2 single crystals and 
derivatives, as well as the ceramic sample. The resistivity of 
the densified ceramic TiS 2 is much higher than the single crys¬ 
tal despite very similar charge carrier concentration. Dashed 
lines show the fits to po + aT^. b) Seebeck coefficient of the 
same samples. 


this sample are in Ref. [H 

Transport measurements were made on single crystal 
platelets of typical dimension 5 x 1 x 0.02 mm. Electri¬ 
cal contacts were made using Dupont 4929 silver paste. 
Electrical resistivity and Hall effect was measured using 
a standard six probe technique. Thermal conductivity 
and Seebeck coefficient were measured on single crystals 
using a custom built steady-state experiment installed 
in a Physical Properties Measurement System cryostat. 
Calibrated heat pipes were used to measure the ther¬ 
mal power entering and leaving the samplei^. Power loss 
(principally from radiation) did not exceed 30% at 300 K. 
This contributes to an additional uncertainty in the ther¬ 
mal conductivity measurements of approximately 10% at 
300 K, which becomes negligible below 150 K. One crystal 
of TiS 2 was also measured using the Harman technique^^ 
and the results are in good agreement. 


III. RESULTS AND DISCUSSION 
A. Electrical Transport 

Electrical resistivity of the samples is shown in 
Fig. [Iji. They show metallic-like behavior with p(300K) 
in between 0.6 and 1 mOcm and resistivity ratio 
p(300K)/p(5K) = 4.5 related to residual resistivity val¬ 
ues at 5K of 0.2 to 0.3mr2cm. The temperature de¬ 
pendence of the curves follows the previously seen p = 
po + aT^ with 1.3 < x < 2 behaviour over a wide temper¬ 
ature range (60-300 K) that has been ascribed to the ap¬ 
proximately equal contributions of inter- and intra-valley 
scattering to the relaxation time^. Values of po and x are 
shown in Table |T] and the fits are included in Fig. [1^. 

A key problem with TiS 2 and derivatives is control¬ 
ling the charge carrier concentration, which arises both 
from deliberate substitution or intercalation, and from 
the tendency for excess Ti^+ to self-intercalate in the 
van der Waals gap. In the following, we establish the 
concentration of substituted and/or intercalated cations, 
niom by reference to well-characterised ceramic series us¬ 
ing the c-axis lattice parameters^J^ii^, and obtain the 
corresponding charge carrier concentration, rie, from Hall 
effect measurements. 

In the case of the TiS 2 single crystal, c = 5.700A 
givesi^ riion = 0.34(8) x lO^^cm”^ and Ug = 1-8(2) x 
10^^cm“^. These values are consistent with a stoichiom¬ 
etry close to Til 025 S 2 , assuming that all of the excess 

Ti4+ 

is intercalated and that all four electrons become 
free carriers. The ceramic sample of the same nomi¬ 
nal stoichiometry presents nearly identical Seebeck coef¬ 
ficient, which we consider to be a more reliable measure 
of the doping level in ceramic materials than rig, as was 
shown in the cuprate superconductorsii. The value of 
|S'(300Ar)| = 250 pV for the single crystals of Ref. d 
corresponds to between 1 and 1.5% excess titanitum ac¬ 
cording to the same ceramic series^^. 

The Coo,o 4 TiS 2 crystal has ng = 1-1(2) x lO^^cm”^ 
and a correspondingly larger |S'(300Ar)| = 106/rV K“^. 
Previous studies of Co intercalation at high concentra¬ 
tions were consistent with the Co^^ oxidation state^^. 
c = 5.683% implies a stoichiometry close to the nom¬ 
inal Coo,o 4 compositioni^, and shows the contraction 
of the lattice characteristic of the strong covalent na¬ 
ture of bonding between the Co^’*' ion and the sul¬ 
phur atoms in adjacent layers!^. By contrast, the pres¬ 
ence of intercalated Ti^+ tends to enlarge the van der 
Waals gap. The measured Ug therefore corresponds to 
nion = 0.55 X lO^^crn”^ and a stoichiometry C 00 . 032 ; 
this is certainly within the margin of error afforded by 
the lattice parameter measurements. There is no sugges¬ 
tion of significant Ti self-intercalation in this case, and 
we speculate that the contraction of the lattice leads to 
a much lower ability of Ti^+ to enter van der Waals gap. 

The highly doped Tao. 1 Tio. 9 S 2 has a much smaller 
|S'(300A')| = 40/iVK“^ and a correspondingly larger 
Ug = 3.7(4) X 10^^cm“^. Ta substitution at the Ti site 
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FIG. 2. a) Thermal conductivity of single crystal composi¬ 
tions based on TiS 2 . The red circles and yellow diamonds rep¬ 
resent the steady state and Harman measurements on nomi¬ 
nally stoichimetric TiS 2 respectively (see text). The electronic 
contribution estimated from the Wiedemann-Franz law (solid 
lines) is subtracted to give the lattice thermal conductivity 
b). The peak height of ^lattice at ~ 80K is reduced as im¬ 
purity scattering increases. The fits to the Debye-Callaway 
model (black solid lines) are dominated by the variation in 
this parameter. The estimated minimum lattice thermal con¬ 
ductivity is also shown. 


provides one free carrier per ion, as can be seen from the 
linear dependence of Ue on doping in Ref. The value of 
|S'(300R')| is close to that obtained for a ceramic sample 
of the same nominal Tao.i composition^, as is c = 5.73324. 
Given this composition, only 1.74 x lO^^cm”^ charge 
carriers are expected. We conclude that the rest of 
the carriers arise from excess Ti intercalation, such that 
the composition is approximately Tao,iTio,g 3 S 2 , giving 
= 1.74 X lO^^cm”^ and = 0.5 x lO^^cm”^. 
The nominal and derived chemical formulae for all of the 
single crystals are summarised in Table HI 


B. Thermal Transport 


The in-plane thermal conductivity (Fig. [2^), partic¬ 
ularly the lattice contribution (Fig. [^b), is a powerful 
probe of riion, especially in the absence of grain bound¬ 
aries. The only previously reported thermal conductiv¬ 
ity value for a slightly self-intercalated TiS 2 crystal^ was 
K(300iT) = 6.8 W with S = -251 /iV cor¬ 

responding to rie = 0.28 X 10^^cm“^. For the latter, 
the electronic part of k was only 0.43 W K“^m“^ against 
Ate/ ~ 1.0 W K“^m“^ for the present TiS 2 and Coo,o4TiS2 
crystals (assuming that the Wiedemann-Franz law is 
valid in this regime for electronic heat transport). Gen¬ 
erally speaking, the thermal conductivity curves follow 
the typical form expected for heat transported mainly 
by phonons: a large quasi-linear rise at low temperature 
as phonons are excited, a peak limited by cation and de¬ 
fect concentration, followed by a continuing decrease at 
higher temperatures as phonon-phonon umklapp scatter¬ 
ing processes become dominant. 

We extract the lattice thermal conductivity (Fig. Oa) 
by subtracting the electrical contribution to the total 
thermal conductivity according to the Wiedemann-Franz 

law, Hilattice = Ktotal “ Kel = Ktotal “ where we 

use the Sommerfeld value of the Lorenz number Lq = 
2.44 X 10“® Kei extracted in this way is usu¬ 

ally an upper bound since a) inelastic scattering pro¬ 
cesses, maximally active at 0 _d/ 5, scatter the thermal 
but not the electrical current effectively and b) if the tem¬ 
perature broadening of the Fermi distribution approaches 
the bandwidth, scattering processes are no longer quasi¬ 
elastic. For the high carrier concentrations in the sam¬ 
ples measured here, the Fermi temperature is in excess of 
1000 K and our measurements remain in the quasi-elastic 
regime^i. 

We analyse niatuce using the Debye-Gallaway model 
with contributions from boundary, impurity, umklapp 
and planar defect scattering: 


cOd/t 


^lattice — 


2'K‘^Vph 


r{x,T) 


4 CL 

xe 


(e^ - ly 


( 1 ) 


where l/r = Vph/d + Auj'^ -I- CuP' and 

X = fiuj/ksT. Vph = 3793 ms“^ is the in-plane sound 
velocity's and d = 1 mm is the relevant limit for bound¬ 
ary scattering, corresponding to the average sample di¬ 
mension. A is the prefactor for scattering by point-like 
defects (lattice imperfections or foreign cations), B and 
b are the parameters describing the umklapp processes, 
and C is related to the density of planar defects. The in¬ 
clusion of C significantly improves the quality of the fits 
at low temperature. The impurity, defect and boundary 
scattering control the shape of the curve at low temper¬ 
atures, while umklapp scattering controls the decline at 
high temperatures. 0 _d = 240 K is used for all samples^!. 
In the case of the ceramic sample, we use d = 33 nm in 
order to compare to the results of Ref. H. 

The results (fits are shown in Fig. [5]3 and parameters 
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Sample 

c(A) 

ne(300K) 

(lO^^cm"®) 

Ilirnp 

(lO^^cm-^) 

po X 

(mflcm) p ~ 

S(300K) 
(pV K”i) 

A 

(A-S-i) 

B 

(A-5s-i) 

b C 

(lO^A-^s-i) 

TiS 2 (crystal) 

5.700 

1.8(2) 

0.43 

0.18 

1.64 

-95 

1200 

0.374 

4.6 

3.61 

Coo.o4TiS2 

5.683 

1.1(2) 

0.55 

0.33 

1.38 

-106 

1850 

0.153 

4.5 

7.21 

Tao,iTio.9S2 

5.733 

3.7(4) 

1.74(Ta), 0.5(Ti) 

0.25 

1.35 

-45 

13400 

0.367 

4.3 

11.6 

TiS 2 (ceramic) 

5.701 

1.1(1) 

1.74 

0.38 

1.71 

-93 

730 

0.455 

2.8 

3.26 


TABLE I. Electronic properties and fit parameters of Debye-Callaway model of the lattice thermal conductivity of TiS 2 and 
derivatives discussed in the text. 


in Table |I| show that the impurity and defect scatter¬ 
ing parameters {A and C) for the single crystals trend as 
expected, increasing as the peak in thermal conductivity 
around 80 K is reduced. The umklapp scattering param¬ 
eters do not show a discernible trend but are all of similar 
magnitude; these are strongly influenced by the high tem¬ 
perature behaviour where uncertainty in the measure¬ 
ments leads to rather large uncertainty in the extracted 
parameters. Additionally, the single crystal data all have 
a similar high-temperature form, with a slight concavity 
near 130 K. This suggests that the variation in niatuce is 
dominated by the changes in the foreign cation concen¬ 
tration, as is supported by the fit parameters. A similar 
trend was observed in a series of Ta 2 Tii_ 2 S 2 ceramics^. 

The A parameter as modelled by Klemens^^ is propor¬ 
tional to three contributions, as shown in Eqn. [21 





( 2 ) 


where nionln is the fractional concentration of substitu- 
tions/intercalants. AM represents the difference in the 
mass of the foreign cation relative to the mass of the re¬ 
placed ion, M. e is a parameter related to stiffness and 
anharmonicity that scales the influence of the size dif¬ 
ference /S.5/5 of the cation relative to the undisturbed 
site. 

Given the proportionality to riion, we can thus calcu¬ 
late the effective A parameter per ion, Aton = Ajuion- 
In the case of two types of cation A = Ai, assuming 
Mathiessen’s rule holds, since different scattering rates 
are additive in EqnjTJ These are shown in Table IIIl Aion 
for Ta'^’*' is approximately 3 times the value for Ti^+ 
intercalants, with Co^"*" intercalants being only slightly 
more effective scatterers than Ti'^+. The parameters 
(AM/M)^ and (AS/S)^ are estimated from the atomic 
masses and observed unit cell volume changes f Table IIIl). 
The very large mass of Ta contributes significantly, but it 
causes a relatively small change in unit cell volume. By 
contrast, Ti and Co intercalants have a smaller mass ef¬ 
fect. For these, we have assumed M = 37.3 atomic mass 
units, the unit cell average, in order to try to make a di¬ 
rect comparison with the Ta case, but this form of Eqn. [21 
is not well suited to intercalants, where a zero mass void 
is replaced by an ion. The volume changes per ion are 
larger for Ti and Co than for Ta, and so the associated 
contributions to the scattering rate are not as small as 
the mass parameters would suggest. 


ion 

Aior^ (lO-^icm^A'-^s-i) 

{AM/Mf 

{AS/Sf 

Ti4-t 

2670 

0.082 

0.31 

Co2+ 

3360 

0.33 

0.15 

Ta^+ 

6900 

8.6 

0.034 


TABLE II. Effective parameters for scattering rate per sub¬ 
stituted/intercalated cation in TiS 2 single crystals. 


The ceramic sample shows a smaller A than the single 
crystal of TiS 2 , but also a much lower Kiattice- The pri¬ 
mary reason for the difference between these two cases 
is the parameter d, which describes scattering by grain 
boundaries. If we take the identical parameters as for 
the single crystal, a value of d = 60 nm produces a very 
similar Kiatuce to the one shown in Fig. [213, showing that 
this is indeed the dominant scattering mechanism. The 
reduction in thermal conductivity on going from crystal 
to ceramic is almost perfectly balanced by the increase in 
resistivity, however, leading to no change in ZT at 300 K. 


IV. CONCLUSIONS 

In conclusion, we have shown using well calibrated con¬ 
centrations of substituted or intercalated ions in single 
crystal materials that the most effective species at reduc¬ 
ing the thermal conductivity in the TiS 2 system is Ta^’*' 
substitution at the Ti site. Ta substitution is particularly 
suitable as it provides only 1 electron per ion, compared 
to 2 and 4 for Co and Ti intercalants, introducing more 
scattering centers for a given charge carrier concentra¬ 
tion. The thermoelectric performance of Ta doped mate¬ 
rials so far studied has however been limited by parasitic 
self-intercalation of Ti. Improvements in thermoelectric 
performance will only be realised once this effect can be 
controlled. Then studies of ceramic TiS 2 and derivatives 
for potential applications can be pursued, as similar per¬ 
formance can be obtained even in the presence of grain 
boundaries. 
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